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Abstract. The spin-statistics connection is derived in a simple manner under the 
postulates that the original and the exchange wave functions are simply added, and 
that the azimuthal phase angle, which defines the orientation of the spin part of 
each single-particle spin-component eigenfunction in the plane normal to the spin- 
quantization axis, is exchanged along with the other parameters. The spin factor 
(— 1)^'' belongs to the exchange wave function when this function is constructed so as 
to get the spinor ambiguity under control. This is achieved by effecting the exchange 
of the azimuthal angle by means of rotations and admitting only rotations in one 
sense. The procedure works in Galilean as well as in Lorentz-invariant quantum 
mechanics. Relativistic quantum field theory is not required. 
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Physics is simple but subtle 
Paul Ehrenfest 

1 Introduction 

The standard method of treating systems of identical particles in present quantum 
mechanics is to require that every wave function or state vector must be either 
symmetric or antisymmetric, that is, multiplied by either +1 or —1 when the labels 
or parameters referring to any two particles are interchanged. There are thus two 
classes of systems, with different collective behaviour of the particles: systems of 
bosons and systems of fermions. These two classes are connected with the spins of 
the particles: all particles which are known to be bosons are empirically found to 
have integral spin, in units of h, while all known fermions have half- integral (i.e., 
half-odd-integral) spin. 

* Extended version of a talk given at the international workshop on theoretical and experimental 
aspects of the spin-statistics connection and related symmetries (SpinStat2008), Trieste, Italy, 21-25 
October 2008. Also published in Foundations of Physics 40 (7), 776 - 792, 793 - 794 (2010) (Special 
Issue on Spin Statistics). This e-print version is the same as the paper- published version except for 
a few places with slight improvements in clarity and precision. 
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Within quantum mechanics the connection with spin could not be derived and 
had to be taken as another postulate. The first derivation was provided by Fierz 
[1] and Pauli [2], who founded it on relativistic quantum field theory. This also 
remained the framework for the papers which in subsequent years refined and gener- 
alized Pauli's proof [3, 4]. Typically, in these papers it is postulated that no negative- 
energy states exist, that the metric in Hilbert space is positive definite, and that 
the fields either commute or anticommute for spacelike separations (locality, micro- 
causality) . Under these conditions it is shown that integral-spin fields cannot satisfy 
the (fermionic) anticommutation relations, and half-integral-spin fields cannot sat- 
isfy the (bosonic) commutation relations. This does not exclude the possibility that 
fields exist which satisfy other commutation relations and show statistics that differ 
from Bose and Fermi statistics. 

In 1965 Feynman in his Lectures on Physics [5, p. 4-3] objected: 

An explanation has been worked out by Pauli from complicated argu- 
ments of quantum field theory and relativity. He has shown that the 
two [spin and statistics] must necessarily go together, but we have not 
been able to find a way of reproducing his arguments on an elementary 
level. It appears to be one of the few places in physics where there is 
a rule which can be stated very simply, but for which no one has found 
a simple and easy explanation. The explanation is deep down in rela- 
tivistic quantum mechanics. This probably means that we do not have 
a complete understanding of the fundamental principle involved. 

The aim of the present paper is to propose such a simple and easy explanation. 

Actually, since 1965 more than a hundred publications appeared deriving the 
spin-statistics connection under different sets of conditions [6]. Reviews are con- 
tained in [7-10]. Many of these publications derive the connection in settings far 
removed from standard (local) relativistic quantum field theory; and they are also 
far from simple and easy. 

Closest to the present approach are those papers that use only quantum me- 
chanics, relativistic or nonrelativistic, and are written in the spirit of Feynman's 
demand for simplicity. These papers nevertheless contain one or several of the 
following restrictions: the wave functions must have special invariance [11, 12], con- 
tinuity [13-17] or symmetry [18] properties, or must lie in special spin-component 
subspaces [19]. The systems considered must be nonrelativistic [13-22], have only 
two spatial dimensions [22], contain only two particles [18, 19], only particles with 
zero spin [13-17] or spin < 1/2 [20, 21], only point particles [23-25], must admit 
antiparticles [26], or the exchange must be considered as physical transportation of 
real objects [11, 12, 23-25, 27]. 

The present proposal is not subject to any of these restrictions. It grew out of an 
attempt to understand the papers by York [28, 29] in the framework of the realist 
interpretation which I developed some time ago [30, 31]. The premises of the present 
proposal are seen when the organization of the paper is considered. 
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In Sect. 2 we start with Feynman's method of superposing transition amplitudes. 
This is shown to be equivalent to symmetrizing or antisymmetrizing only the final 
but not the initial wave function in the transition amplitudes. Moreover, in line 
with another proposal by Fcynman, symmetrizing or antisymmetrizing the final 
wave function, which for two particles means addition or subtraction of the original 
and the exchange wave function, is replaced by the "addition postulate", which 
admits only addition. 

The minus sign in the superposition of fermionic wave functions arises from the 
construction of the exchange function. In Sect. 3 we point out that it is important to 
consider the azimuthal spin angle which defines the orientation of the spin part of 
a single-particle spin-component eigenfunction. The angle x is also exchanged, but 
requires a special treatment because it expresses the well known spinor ambiguity 
when the spin component m is half-integral. This means that we cannot know which 
of the two possible values of the function with the exchanged x has to be chosen. 

In Sect. 4 it is shown that the ambiguity can be overcome by effecting the ex- 
change of X by way of rotations and by admitting only rotations in one direction, 
either clockwise or counterclockwise. In Sect. 5 the rotations leading from the orig- 
inal to the exchange function are explicitly carried out, and it is shown that the 
exchange function thereby acquires the desired spin factor (—1)^*. It is thus the ro- 
tation group, a subgroup of both the Galilei and the Lorentz group, that determines 
the type of statistics. In standard quantum mechanics it is the permutation group 
that does this: its one-dimensional representations are associated with Bosc and 
Fermi statistics, and its other representations with "parastatistics" . The rotation 
group in our approach leads only to Bose and Fermi statistics, and there is no reason 
to suggest experiments in search of particles with parastatistics. 

In Sect. 6 the proof is extended to N particles, and in Sect. 7 to particles with 
different spin components. Finally, in Sect. 8, using the properties of helicity func- 
tions, it is pointed out that the proof also holds in relativistic quantum mechanics. 

2 Adding up transition amplitudes 

It is remarkable that wave-function symmetrization (or antisymmetrization) is never 
mentioned in the Feynman Lectures on Physics [5] and yet the same physical sit- 
uations as in the other textbooks, which do use symmetrization, are acounted for. 

Fcynman basically considers transition amplitudes: when two transitions cannot be 
distinguished in principle from each other, the amplitudes, rather than the proba- 
bilities, have to be superposed [5, pp. 1-10, 3-7,. ..,4-3]. The superposition includes 
the phase factor exp(i(5) 

f = f{e) + e''f(j,-e) (1) 

in Feynman's notation. In line with standard quantum mechanics the phase factor 
in the Lectures is eventually postulated to be -|-1 or — 1. 

What does this mean in terms of (Schrodinger) wave functions? The transition 
amplitude for the transition from ^'a(l, 2) to 2) is given by the scalar product 

/= (%(l,2),M/a(l,2)). 
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These scalar products are basic elements of quantum mechanics because any expres- 
sion of physical significance, that is, an expression which yields the probability of 
a result of a measurement, is formulated in terms of absolute squares of transition 
amplitudes 

|/P = |(*b,*a)P. 

Now, in systems of identical particles the transition from '^^,{1, 2) to ^'b(l, 2) cannot 
be distinguished in principle from a transition from ^'a(l,2) to '^\y(2,l), where 
'J'l-,(2, 1) denotes the exchange function, in which the labels or parameters referring 
to particles 1 and 2 have been exchanged. Thus, in terms of wave functions (1) takes 
on the form 

/ = {^^{1, 2), ^'a(l, 2)) + e''(^^{2, 1), *a(l, 2)) 

= (*b(l,2)+e'^*b(2,l),^a(l,2)), (2) 

and with e''^ = ±1 

/=(*b(l,2)±*b(2,l),*a(l,2)). (3) 

Here only one of the two functions in the scalar product is a superposition of the 
original and the exchange function, and no normalization factors appear. This is 
the form that is exclusively used by Feynman. 

This is to be compared with the standard method in quantum mechanics, where 
the amplitude in the above case is written as 

/= (^[^b(l,2)±*b(2,l)],^[^a(l,2)±*a(2,l)]). (4) 

Here both functions in the scalar product are superpositions of an original and an 
exchange function, and both functions are normalized, by means of the factors 1/V2. 

Feynman's formula (3) and the standard formula (4) are equivalent provided the 
phase factor exp(i(5) is restricted to the values ±1. The proof of this is based on 
properties of the symmetrizer S and the antisymmetrizer A, 

S = (1/iV!) , A= (1/iV!) E 

where N is the number of identical particles {N > 2 here included), Pq, is a per- 
mutation, and ka is the number of transpositions (interchanges) that make up the 
permutation Pa- The operators 5" and A satisfy the relations 

S^ = S = S'^ and A^ = A = A'^, (5) 

the proof of which can be found in [32, pp. 1383 - 1385]. 

With these relations the standard formula (4) in the Fermi case, when written 
for AT-particle functions 



/ = (VM^^'b(l, 2, . . . , iV), VM^*a(l, 2, . . . , iV)) 



(6) 
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becomes 

f = m (A^A^h, *a) = N\ (A^b, *a) = ( ^^(-l)'="Pa*6, ^a) , (7) 

which is Feynman's formula (2) with e"^ = — 1, when written for A''-particle functions. 
In the Bosc case the equivalence is obtained in an analogous way. 

When the factors in front of the P^s in S* and A, that is 1 and ( — 1)*^", respec- 
tively, are replaced by more general factors inspection of the proof of relations (5) 
in [32] shows that the proof no longer goes through, and with this the equivalence 
between the standard formulas (4), (6) and the Feynman formulas (2), (3), (7) no 
longer holds. Such a more general factor appears however in the course of our treat- 
ment when it comes to one-particle wave functions with different spin components 
{ria, formula (28)). We thus have to make a choice. It is Feynman's formulas that 
we choose because we found that in contrast to the standard method the Feynman 
method allows us to incorporate a derivation of the spin-statistics connection and 
yet finally to arrive at the approved expressions of physical significance. 

Physically, Feynman's method recommends itself because it is really only nec- 
essary to consider physically significant expressions, and the scalar products are 
closer to these (|/P) than are the wave functions (cf. e.g. [33, 34]). In any case 
it appears reasonable to number the particles in the first function arbitrarily and, 
assuming that the numbering gets lost in the transition, to add up all possibilities of 
numbering in the second function. In (anti)symmetrizing only the second function 
the Feynman method is less restrictive than the standard one. This may be the 
deeper reason why the Feynman method allows for a derivation of the spin-statistics 
connection, while the standard method does not. 

In choosing the Feynman method we shall see that finally in all expressions of 
physical significance the more general factor rja reduces to (— i)^**:^^ This allows us 
to write e"' = (—1)^'' in Feynman's transition amplitudes (2) and (3), and (—1)^" = 
(^_l^2ska ]^jg amplitude (7). These factors are restricted to the values ±1, like 
the factors in S and A. With this the equivalence between the standard and the 
Feynman method will be recovered, with the additional bonus that the Feynman 
method has enabled us to establish the connection with spin. 

In his own attempt to derive the spin-statistics connection Feynman [27] more- 
over suggests that we may 

take the view that the Bose rule is obvious from some kind of under- 
standing that the amplitude [s] in quantum mechanics that correspond 
to alternatives must be added. 

That is, he is proposing to start a priori with exp(i(5) = everywhere. We follow 
also this proposal by Feynman and just add up the transition amplitudes. That is, 
we start from 

/= (*b(l,2) + ^'b(2,l),^'a(l,2)) 
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in the case of two particles. In the case of N particles this obviously generalizes to 
/ = (Ea ^«*b(l, 2, . . . , iV), *a(l, 2, . . . , iV)) , 

where is a permutation of the parameters referring to the single particles. The 
sum extends over the N\ possible permutations, including the identity I. The sum 
J2a -fa^b = ^bS' extension of the sum of the original and the exchange func- 

tion from two-particle to AT-particle systems and corresponds to the symmetrized 
function in standard quantum mechanics. 

3 The special parameter 

In order to present the essential points in a simple way we begin by considering a non- 

rclativistic (Galilean) system of two identical particles of spin s {S'^ ip = s{s + 1) ip) 
described by a (Schrodinger) wave function which is a product of two normalized 
one-particle wave functions 

* = tp^^\a, m) tp^'^\b, m) , 

and the one-particle wave functions are eigenfunctions of the operator of the spin 
component with respect to an arbitrary but common spin-quantization axis. More- 
over, both functions belong to the same eigenvalue m. These restrictions will be 
removed in Sects. 6 to 8. 

The single-particle wave functions are functions of the variables x, y, z, t. The 
functions arc determined by the parameters a, b and m, where a and b stand for the 
sets of parameters that, together with m and x (below), allow for a complete account 
of all aspects and degrees of freedom of the single-particle systems. Such a set of 
parameters includes e.g. mass, charge, total spin, centre, expansion coefficients etc. 
Mass, charge, and total spin are of course the same for identical particles and their 
exchange has no effect. An alternative notation would be tp{a,m,x^^\y^^\ z^^\t), 
where the labels in parentheses, (1) and (2), distinguish the particles in the formal- 
ism. We have suppressed here the variables and have put the particle labels directly 
at the function symbols. In Sects. 2 to 6 the eigenvalue m is always the same and 
will also be omitted from the notation. Thus we write 

= ^(i)(a) . (8) 

Among the parameters of the wave function there is one that requires special treat- 
ment in the construction of the exchange function because it may lead to double- 
valued wave functions. The reason is that the spin parts of the wave functions, 
while belonging to one and the same m, may still differ from one another by a rota- 
tion about the spin-quantization axis. In other words, each spin part has a definite 
orientation in a plane normal to the common spin-quantization axis, defined by an 
azimuthal angle x, counted from some arbitrary reference direction. A complete 
spin-quantization fram,e rather than only a spin-quantization axis is involved. The 
angle x is kept out of the set a (and b) and is exhibited explicitly. Each function 
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can have its own angle, but the particular values do not matter. The values of x are 
restricted to the interval [0, 2tt]. 

The specific form of the parametric dependence of the spin-component eigen- 
function on x is given by the factor 

exp(imx), (9) 

so that 

Xa) = exp(imxa) (10) 

The angle x appears "only" in a phase factor, and when this factor is an overall 
(global) phase factor it is without any physical significance and can be omitted. 
However, it must be taken seriously if it is to become part of a superposition, thereby 
determining a relative phase and thus becoming physically significant [32, pp. 219, 
220]. This is what happens in the present approach: the angle x is exchanged, 
though in a specific way, along with the other parameters, and in Sect. 5, when 
the original and the exchange function are superposed, x becomes instrumental in 
determining the relative phase between these functions. 

The exponential factor (9) expresses the spinor ambiguity: in the case of half- 
integral m it is -1-1 for X = and — 1 for x = 27r. The sign change under a full 
rotation holds for every rotation axis. The factor (9) refers to the special case of 
rotations about the spin-quantization axis. This docs not lead out of the subspace 
of functions with the same m. All this is standard quantum mechanics [32, pp. 694, 
703, 985, 986]. 

Note that the rotation of the orbital angular momentum part of a wave function 

can be expressed by changing the value of an already present spatial variable {ip — t- 
if + f' [32, pp. 681, 682]), but that an additional parameter (x) is needed to express 
the behaviour of the spin part of a wave function under a rotation. 

The letter x, rather than the customary ip, is used in order to emphasize that it 
is the spin part, not the orbital part, which is concerned, that the spin-quantization 
axis need not coincide with the z-axis, and that the angle x is not a variable of 
the one-particle wave function, as r and if are. Rather, the dependence on x is 
a parametric dependence, like that on m and the other parameters in a and b. 
Therefore the application of a differential operator like —ihd/dx, analogous to the 
z-component of orbital angular momentum, does not make sense for the spin part 
of the wave function. 

Thus, with the explicit appearance of x (8) becomes 

^ = ^^%,Xa)^^^\b,Xb). (11) 

In constructing the exchange wave function in traditional quantum mechanics it is 
irrelevant whether we exchange the particle labels (1), (2) or the function param- 
eters a, b. In our construction it is no longer irrelevant, and it is the exchange of 
the parameters that must be chosen. Thus we replace a by 6 and vice versa in the 
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original wave function (11). But because the parameter x expresses the spinor ambi- 
guity the exchange of Xa with Xh aiid vice versa cannot be done in such a simple way. 

4 Controlling the spinor ambiguity 

The special feature with the factor (9) is that x is an angle, so that we may go from 
some particular value Xa to some other value Xb in two ways, either clockwise or 
counterclockwise. In the case of half-integral m one way leads to a different wave 
function at Xh than the other, the two functions having different signs. In other 
words, the value of the function at Xb then depends not only on the value of Xb but 
also on the path leading from Xa to Xb- This leads to double- valued functions and 
represents another aspect of the spinor ambiguity. 

One may imagine the function exp(imx) with half-integral m to lie on the two- 
sheeted Riemannian surface of the function yfz [35, 36], where one sheet carries only 
one set of function values. The clockwise path from Xa to Xb always ends up in a 
different sheet than the counterclockwise path. Or one may imagine a Mobius band, 
where on the first round trip over the band one set of function values is met, and 
the corresponding other set on the second round trip. In fact, devices like twisted 
ribbon belts [27, p. 58], contortions of an arm holding a cup [27, p. 30] and others 
[37-42] are similar to the Riemannian surface and the Mobius band in that they 
construct an indicator of whether we are in the first or in the second turn, and 
in that they return to the original situation after the second turn. For integral m 
(including s = m = 0) the Riemannian surface has only one sheet and no ambiguity 
arises. 

Now, when adding the original and exchange wave functions the functions must 
be uniquely defined. This is not the same as the general requirement that wave 
functions be single- valued. Single-valucdness can only be required for measurable 
quantities such as transition probabilities or expectation values, but not for the wave 
functions themselves [43] . In many textbooks it is nevertheless invoked for the wave 
functions themselves, in particular for justifying the restriction to integral values of 
m for orbital angular momentum. The real justification of integral m here rests on 
group representations and properties of observables [44]. 

Our case is different because we are concerned with the procedure of construct- 
ing one wave function by superposition of others, formally similar to interference. 
The demand for removing the spinor ambiguity is in line with the demand for re- 
moving the ambiguity known as exchange degeneracy, that is, to the fixation of the 
coefficients in the superpositions in Sect. 2. 

Now, according to what has been said above the spinor ambiguity is removed 
(i.e., kept under control) if we make a choice between the two possible paths from 
Xa to Xb-, that is, if we exchange the xs by way of rotations and decide to make all 
rotations in one sense only, either clockwise or counterclockwise. 

In the language of group theory the clockwise and the counterclockwise way from 
Xa to Xb correspond to paths of different homotopy classes (e.g. [45]). So our choice 
means that we are admitting only paths of the same homotopy class. 
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5 Constructing the exchange function 

We are now ready to take the decisive step. We want to construct the exchange 
function from the original function (11), not by simply replacing Xa by Xb iii the 
wave function -^^^^ and Xb by Xa in the wave function ip^"^^ (as is done with the other 
parameters, a and b), but by continuously rotating the spin part of the functions 
from Xa to Xb and from Xb to Xa respectively, with due consideration being given to 
the paths connecting Xa and Xb- 

Thus, we start from formula (11) where the as and 6s have already been ex- 
changed, but the xs have not: 

ij^'\b,Xa)^P^%,Xb). (12) 

We then rotate the function tp^^Xb, Xa) from Xa to Xb- We take the counterclockwise 
sense of the rotations, and we assume Xa < Xb and m > 0. In order to get from 
Xa to Xb we then have to run through Xb — Xa- This yields the rotation factor 
exp(im(x6 — Xa)) and we obtain 

i;^'\b,xb) = e'"'^^''-^'^^^'\b,xa)- (13) 

Likewise, rotating the function ip^'^Xa, Xb) counterclockwise from Xb to Xa means 
that wc have to run through 2tt — (xj^ — Xa)- This yields the rotation factor 
exp(im(27r + Xa — Xb)) and we obtain 

V'(2)(a,x„) = e'-(2-+x.-x.)^(2)(^^^^). (14) 
Inserting (13) and (14) into (12) then yields the exchange function 

Fx^^'\b,Xb)^^'\a,Xa) (15) 

with 

F = Q-^"^iXb-Xa) Q-im{27T+Xa-Xb) — g-im27r _ (_i)2m _ (_i)2s^ (-|^g) 

where for the last equality we have used the fact that s and m are either both 
integral or both half-integral. Had we chosen the clockwise sense we would have 
obtained F = exp(-|-im27r), which is also equal to (—1)^*. The same result obtains 
when m < or when Xa > Xb- The case Xa = Xb is of statistical weight zero and 
can be neglected. 

Finally, adding the original function (11) and the exchange function (15) we 
arrive at 

*S = i^^%,Xa) i^^''\b,xb) + (-1)'^ i'^'\b,Xb) i^^^\a,Xa)- (17) 

^'S need not be normalized in Feynman's method, as emphasized after formula (3). 
The angle x and the rotations become effective only in the procedure of exchang- 
ing the parameters of the wave functions. In this procedure the angles Xa and 
Xb in the original and the exchange wave function are related in such a way that, 
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although Xa and Xb may be randomly distributed in the original function, the resul- 
tant factor, (—1)^'', is independent of Xa and Xb- Thus, once exchange and addition 
are accomplished, we may put the angles x back into the parameter sets a and 
b. Or we may, according to formula (10), write the functions in (17) in the form 
ip^^\a, Xa) = exp(imxa) tp^^\a) etc. Either term in (17) thereby receives the same 
factor 

exp(im[xa + Xb]), (18) 

which can even be put as an overall phase factor in front of the sum, where it can 
be omitted, that is, absorbed in the general arbitrary phase factor connected with 
every wave function. 

Thus we are returning to the standard form of the wave functions, which do not 
explicitly exhibit the dependence on x- 

*S = tp^%) il^^'^Xb) + {-if il;^^\b) ijj^'^Xa). (19) 

There is some formal analogy with interference between two parts of a split wave. 
One part is left unmodified [wave function (11)], the other is subject to a phase shift 
[exchange, wave function (15)], and then the two are recombined [wave function (17) 
or (19)]. 

When (19) is inserted into Feynamn's formula (2) with ^'(,(1, 2) -|- e^^ *fe(2, 1) =: 
we see that e'^ = (—1)^^, and with this we have reached our goal for the con- 
sidered class of functions: the factor (—1)^'' is no longer postulated but is derived 
in a simple way from basic principles. This factor yields +1 (bosons) for integral 
s and —1 (fermions) for half-integral s, and this is the desired connection between 
spin and statistics. 

6 General case. Equal spin components 

We begin now to remove the restrictions imposed on the wave function in the pre- 
vious sections. In the present section we first remove the restriction to two particles 
and then to functions of product form. The A/"-particle functions of product form 
are 

*b = 'll^'-^Kur^ , Xh ) V'^^^Wra , Xta ) " " " Vriv i Xt^ ) 

with instead of a,b. The tjj^^\urf,Xti) all belong still to the same m, 

which is therefore dropped from the notation. The symmetrized function is 

*bS' = V'^'^^n , Xt J • • • i^^^'Kur^ ,Xt^ )- (20) 

The index S' (with prime) is to indicate that Pa in (20) permutes the parameter 
sets among the one-particle functions but does not permute the angles Xu ■ The 
permutation of the angles will be effected separately, by way of rotations. As any 
permutation can be written as a product of a number of transpositions, the term 
Pa '4'^^\uri , Xti ) ■ ■ ■ '^^^\uri^ , XtN ) differs from the term with P„ = by a number kce 
of transpositions. When the x rotations are applied, as described in the preceding 
sections for the case of two particles, every single transposition yields the factor 



- 11 - 



F = (—1)^* in front of the term with interchanged parameters, independent of the 
angles X- Hence transpositions yield the factor (_i)2sfca^ The function (20) then 
changes into the superposition function (symmetric or antisymmetric) 

*bS = Ea(-l)'''"^« V'^Vi,Xti)---^^^W,Xt.)- (21) 

The index S (without prime) instead of S' (with prime, as in (20)) is to indicate that 
the exchange of the angles xu by means of rotation is included. In (21) therefore 
permutes the pairs {uri,Xti}- The single functions (— l)^*''^"Po, , Xti ) " " ' 

ip^^\ur^,XtN) Pa ¥^ I are the extensions of the exchange function from two to 
A'" particles. 

As in the case of two particles, the one-particle functions in (21) may now be 
written as ip^'^\uri,Xti) = exp(imxtj ?/'''*\^r-i)- This yields the same factor 

exp(imEfcLiXt;,) 

in front of every permutation operator in (21). This factor can thus be drawn 
out of the permutation sum and becomes again an overall phase factor, which can 
be omitted. We can thus drop the x^. s from the notation in (21) and simply write 

*bS = J^.i-l f^'-Pa ^^\ur,) ■ ■ ■ V^n^r.). (22) 

Now, if s is an integer, then (— i)2«*:a _ foj- a,ny ka, and 'fbs is symmetric 
(bosonic). If s is a half-integer, then (— l)2sfca — _\ fgj. qjJ^j /^^^ g^j^^j j^-y ^qt^ even 
fea, and \I'i3s is antisymmetric (fermionic). And this holds for each one of the 2s+l 
values of m. 

Now we remove the restriction to wave functions of product form, but still with 
equal ms. The general iV-particle function then is 

V'^^W^Xti) • • • V'^^Wjv!Xtjv)> (23) 

r-i,...,r-Ar,ti,...,tjv 

where the sum (or integral) over the rs and ts goes over a possibly infinite number 
of values. Permuting the parameter sets Ur^ among the one-particle functions and 
permuting the angles xu by way of rotations now results in 

^bs = E ^,....t,..t.E(-i)'''"^«^^Vi,xti)---V'^'^W,xt.)- 

ri,...,rjv,ti,...,tAr a 

(24) 

The sums '}2a be considered as basis functions in the subspace of wave functions 
with equal ms. They are symmetric or antisymmetric. The antisymmetric ones are 
Slater determinants, on which the Pauli exclusion principle is based. 

From the consideration in the preceding section (cf. (18)) it follows that the 
permutation sum in (24) can be written as 

exp(im[x*i + • • • + XtA) Ej-l)'''"^a ^^'W) • • • V^^'^W)- (25) 
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So 

*bS = Wi-TNti-tN exp (im[xti + • • • + XtA) 

ri,...,rN,ti,...,tM 

X x:(-i)2^^"PaV'(Vi)---V'^^W)- (26) 

a 

The exponentials can still be drawn out of the permutation sum, as in (17) and (21). 
They can no longer be drawn out of the sum over the rs and is in (24) and cannot 
become an overall phase factor. This is not essential. At this stage the identity of 
the particles has already been duly taken into account, and the sum over the rs and 
ts has no longer anything to do with it. The appearance of the exponentials under 
the sum only means that the angles x a-re explicitly exhibited, rather than being put 
into the sets (cf. the remarks after Eq. (17)). 

The connection between spin and statistics is thus proved for general nonrela- 
tivistic iV-particle functions composed of one-particle functions with the same spin 
component. 

7 Diflferent spin components 

When admitting one-particle wave functions which belong to different spin compo- 
nents m, the factor F in (16) of Sect. 5 has to be replaced by a different factor, F^. 
This factor is obtained when the transposition procedure of Sect. 5 between (12) 
and (16) is repeated with (12) being replaced by tp^^\b,mb,Xa) '4^'^'^\o-,fnaiXh)- The 
result is 

F^ = i-lf' exp(-i(m„ - rribXxa - Xb))- (27) 

As any permutation can be written as a product of a number ka of transpositions, 
the factor (— l)2sfca ^21) to (26) is then replaced by r)a, which is a product of F^s: 

rja = i-lf''- n^"^ exp(-i(m„ - m^Xxa - Xb)il - 5fc„o) ) • (28) 

The values that the parameters iria, Xa, t^b-, Xb assume in the various factors of the 
product are those of the particular one-particle functions on which the transpositions 
in Pa operate. The factor (1 — 5fc^o) with the Kronecker delta is to make the product 
1 when ka = 0, that is, when Pa = I- As we shall see, the exact form of rja does 
not matter. What matters is first that in the case of equal ms the factor Fy. reduces 
to F, and rja to (— gQ ^j^g^^ g^.^ back to Sects. 2 to 6, and second that 
in any case \ria\ = 1- What we have to show is that even in the case of different 
ms the factor r)a is effectively (— l)2sfea ^^^j ^j^g angles x will disappear from the 
permutation sum. 

The general A''-particle wave functions now are (cf. (23)) 

Ml,---,iV) = 

= J2 r•l,...,r^^ a^i...rjv^i-3jvti-tjv V'^^ W , , Xti ) " " " V'^^ , r"s;v ' Xtjv) (29) 

SI, ... , SiV 
tl, ■ ■ ■ ,tN 
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$b(i,...,iv) = 



(Tl , . . . , (TJV 

n, • • • , Tat 



where the sum over the Sj s and aj s goes over the 2s + 1 possible values of the spin 
component. We then have 

a 

= bp^...r^YriaPai'^\up^,ma^,XT,)- ■ ■i'^^\upf,,mai^,XTM)- (31) 

pi,...,TN a 

With a view to an expression of physical significance we consider the transition 
amplitude 

/= ($bS,^a) =E^Pi-r^E«n-*;vrMpar), (32) 

where 

a 

)), (33) 

and the Pa permute only the sets {pi,ai,Ti}, not the sets {ri,Si,ti}, among the 
one-particle functions. The point is that the one-particle functions with different 
spin components are mutually orthogonal, irrespective of the us and ^s- Thus, the 
scalar products between these functions are all proportional to Kronecker deltas, 
and the terms (33) become 

'P{rstpaT) ~ VaPoi'^o-isi ' ' ' I^omSn ^<J\s\ ' ' ' ^omSni (3^) 
a 

where k^.^, = "i<7,, XrJ, ip^'\ur^,msi,XtS) ■ 

The only non-zero terms in the sum are those where in each Kronecker delta 
the pair of indices consists of equal numbers, ai=Si for each i, although different 
Kroneckers may have different pairs. 
There are three types of these terms: 

(i) Terms where all fjjS, and hence all SjS, are equal. These cases are those 
already solved in the preceding sections, yielding the desired result rja = (— l)^*^". 

(ii) Terms where all (jjS, and hence all SjS, are different, though the {si} and the 
{ai} are from the same set of numbers {N < 2s -|- 1). In these terms the sum 
over the permutations reduces to one single member, in which the permutation Pa 
of the sets {upi,Tna^,XTi} achieves that ai=Si for each i. The term T then reduces 
to 

'P{rstp(TT) ^ii 

• "0^ \upj^,msj^-iXTN^ 1 
i;^^\ur„ms„Xh) ■ ■ ■ i^^^\urj,,ms^,Xtj) = r?*^ (^-b, ^-a), (35) 
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where *a and ^'^ are product functions like those in (29) and (30) but without the 
sums, and in the special case o"j=Sj. In the transition amplitude (32), (33) the rj^^s 
of the terms (35) are overall factors in front of the permutation sum {T^^rstpar)) and 
only contribute to the sum over the p,a,r,r,s,t in (32), which is not essential, as 
explained after formula (26). 

In the special subcase that the product functions 'I'a and already repre- 
sent the total wave functions the transition amplitude (32) is equal to (35), that 
is, to one term only. The same then holds for the probability of the transition, 
|(*l3,*a)p5 ilViil = 1)) and there are no interference terms involving different tran- 
sitions. This is as in a system of distinguishable particles. We have here a general- 
ization of a a well known result for two particles [5, p. 3-12], [32, pp. 1407, 1408], 
[46]. 

(iii) Terms which consist of sets with both equal and unequal crjS. In these 
terms the sum J2a ^^^^ the permutations of the sets {wp . , mg-i , Xn } reduces to one 
member in which the permutation Pa achieves that ai=Si for each i, plus members 
in which Pa is supplemented by transpositions between the sets with equal ajS. Pa 
is associated with the factor 7]*^, which depends on the angles x of the transposed 
sets with different CTjS, comparable to the Pa in point (ii). The factor ry*^ is the same 
in all the surviving members. It can thus be drawn out of the permutation sum J2a 
and becomes unessential. 

Thus in all expressions of physical significance the general factor rja of (28) can 
be replaced by the factor (— l)2sfea^ in which there is no longer any dependence on 
X and which can only assume the values ±1. Thus finally we see that indeed the 
standard and the Feynman method are equivalent, as announced in Sect. 2. 

Another point deserves to be mentioned: the term Tf^rstpar) of (34) is zero if the 
values of the {sj} are from a different set of numbers than those of the {cxj} because 
in this case no permutation can achieve that ai=Si for all i. The terms (34) thus 
seem to be restricted to spin-independent interactions. Such interactions could how- 
ever be included by considering the possible wave functions that result from such 
an interaction as intermediate functions between ^'a and ^bS; and by forming 
the sum over the product of the amplitudes X^i(^bS' ^l){^b *a) in the case that 
the (situations associated with the) are not observed, or over the product of the 
respective probabilities in the case that the are observed. 

8 The relativistic domain 

The derivation of the spin-statistics connection presented so far evidently does not re- 
quire relativity theory. Can it be extended into the relativistic domain? In Lorentz- 
invariant theory spin and orbital angular momentum are no longer separately con- 
served quantities, and the two are in general mixed up in a complicated way. There 
are however functions which are eigenfunctions of the spin-component operator only, 
with no admixture of orbital angular momentum: the helicity functions [47]. A he- 
licity function describes a free particle with definite non-zero linear momentum and 
is an eigenfunction of the operator of the spin component with respect to an axis 
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that is parallel or antiparallel to the direction of the particle's momentum. Thus 
we may replace the previously discussed eigenfunctions of the operator of the spin 
component along a fixed direction by the helicity functions. Helicities are invariant 
under ordinary rotations (involving spin and orbital part), and the rotation opera- 
tors commute with the permutation operators, so we may express the momentum 
eigenfunctions which have their momenta in arbitrary directions by suitably rotated 
eigenfunctions with momenta in one common direction (cf. [47, pp. 407, 408]). For 
these functions we can define a common reference direction for the angles Xi 
then construct and add up the functions with the permuted parameters in the pre- 
viously described way. This works not only for momentum eigenstates, i.e. plane 
waves, but also for linear superpositions of plane waves, i.e. wave packets. 
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